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Abstract 

Context. We found in previous studies that standard Simple Stellar Population (SSP) models are unable to describe or explain the 
colours of Galactic open clusters both in the visible and in the NIR spectral range. The reason for this disagreement is the continuous 
nature of the stellar IMF in clusters which is the underlying assumption in the SSP models. In reality, the Galactic open clusters are 
scarcely populated with the brightest stars responsible for integrated fluxes. 

Aims. In this study, we aim at constructing discrete SSP-models which are able to adequately describe the observed magnitude-, 
colour-, and mass-to-luminosity-ratio-age relations of open clusters by including a number of rarely considered effects. 
Methods. We construct a numerical SSP-model, with an underlying Salpeter IMF, valid within an upper m u and lower mi stellar 
mass range, and with total masses M c = 10 2 . . . 10 4 itiq typical of open clusters. We assume that the mass loss from a cluster is 
provided by mass loss from evolved stars and by the dynamical evaporation of low-mass members due to two-body relaxation. The 
data for the latter process were scaled to the models from high-resolution N-body calculations. We also investigate how a change 
of the m; -limit influences magnitudes and colours of clusters of a given mass and derive a necessary condition for a luminosity and 
colour flash. 

Results. The discreteness of the IMF leads to bursts in magnitude and colour of model clusters at moments when red supergiants 
or giants appear and then die. The amplitude of the burst depends on the cluster mass and on the spectral range; it is strongly 
increased in the NIR compared to optical passbands. In the discrete case variations of the parameter mi are able to substantially 
change the magnitude-age and M/L-age relations. For the colours, the lowering of mi considerably amplifies the discreteness effect. 
The influence of dynamical mass loss on colour and magnitude is weak, although it provides a change of the slopes of the considered 
relations, improving their agreement with observations. For the Galactic open clusters we determined luminosity and tidal mass 
independent of each other. The derived mass-to-luminosity ratio shows, on average, an increase with cluster age in the optical, but 
gradually declines with age in the NIR. The observed flash statistics can be used to constrain mi in open clusters. 

Key words. Galaxy: open clusters and associations: general - Galaxy: stellar content - Galaxies: fundamental parameters - Galaxies: 
photometry - Galaxies: starburst - Galaxies: star clusters 



!>• ■ 1. Introduction our pre vious work. Thes e are cluster ages (see Kharchenko et alJ 

!L> _ , , , 2005a! b), tidal masses (Piskunov et al. l l2007l). and integrated 

; Star clusters are the best tracers of star formation and evolu- photometric magnitud es (Kharchenko et al. 2009, referred here- 

tion of Galactic populations, both in the Milky Way and in ex- a ^ ter as p a „ er 
ternal galaxies. In their majority, extragalactic clusters are ob- 
served as unresolved objects. To determine their parameters one The comparison of the observed cluster colours with 

frequently uses Simple Stellar Population (SSP) models based the prediction of the standardSSP-models GALEV 

on methods of evolutionary synthesis. As a benchmark of such (Anders & Fntze-v. Alvensleberj |2003|) and Starburst99 

models, it is reasonable to use the data on resolved clusters d Vazquez & Leithererj | 2005|) revealed significant discrepancies 

where the physical parameters can be determined independent for clusters y° un 8 er than 1 °y r - The observed clusters usually 
of the integrated light. Such clusters should be resolved and ob- bluer than the models P redlct ' but from time to time they 

servable down to faint magnitudes. Obviously, the best candi- &w a sudden "red colour flash". Our analysis ( |Piskunovet al. | 

dates for this kind of investigations are open clusters in nearby l 2 ^, Paper II) has shown that these discrepancies occur when 

galaxies or in the Milky Way itself tne discrete nature of the IMF is neglected in the models. 

We apply this approach to the local population of Galactic In this paper we consider the properties of SSP- 

open clusters. The basic data on open clusters are taken from models and try to achieve compatibility of the models 

with open cluster data by applying the well known ap- 

Send offprint requests to: R.-D. Scholz proach of stochastic modelling of star clusters (see e.g. 
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Figure 1. The empirical colour-magnitude diagram based on 
data from the ASCC-2.5 of the most probable cluster members 
used for the classification of 650 open clusters. The solid b lack 
curve is the empirical ZAMS of Schm idt-Kaler et all (Il982l) . the 
solid cyan curve is the TAMS taken from the Padova models, 
and the broken (magenta) curve separates the Red Giant do- 
main from post-MS stars. The members located to the left of 
the TAMS are classified as Main Sequence stars (marked in 
blue colour), those located to the right of the broken line as 
Red Giants (red), and stars in between (green) are cluster mem- 
bers which just have left the Main Sequence but are not yet red 
enough (Post-MS stars). 



Bruzual 2002 ) adopted in many cluster s t udies (see e.g 



Girar di et al.lll995b ICervino & Luridianall2004t iRaimondo et al.l 
2005 i lPopescu & Hanson! 2010t Fouesneau & Lanconl2010l) . In 
order to make our models as close to reality as possible, we in- 
clude a number of effects which might be important in the evo- 
lution of open clusters of low or moderate masses. The paper 
analyses the evolution of luminosity, colour, and the mass-to- 
luminosity ratio with cluster age, comparing SSP models based 
on a continuous IMF with models based on a discrete IMF. 

The paper has the following structure. In Sect. [2] we briefly 
describe the data on open clusters and adopted definitions. 
Sect. [3] gives the recipes for the construction of realistic SSP- 
models. In Sect.|U|5] and[6]we discuss the evolution of the SSP- 
models in magnitude-age, colour-age, and M / L-age diagrams 
and compare the model results with observations. In Sect. [7] we 
derive the condition for producing a luminosity flash in a star 
cluster and use this condition to derive additional properties of 
Galactic open clusters. In Sect.|8]we summarise the basic results 
of this study. 

2. Data and definitions 

In this study, we consider a sample of 650 local open clus- 
ters identified in the ASCC-2.53. In the latest version of 
the ASCC-2.5 catalogue, NIR magnitudes have been included 
which were taken fro m the 2MASS Point Source Catalogue 
dSkrutskie et alj 120061) . For each cluster, the basic astrophysi- 
cal parameters (age, distance etc.) were determined from ho- 
mogeneous photometric and astrometric data by taking into 
account the membership probabilities of cluster stars (see 
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Figure 2. The empirical colour magnitude diagrams of 650 open 
clusters identified in the ASCC-2.5 in the optical (left panel) 
and the NIR (right panel) photometric passbands. The symbols 
follow the classification scheme shown in FigQ] The blue cir- 
cles mark the clusters with Main Sequence stars only, the green 
crosses are those clusters, where Post-MS stars are also found, 
and the red pluses indicate clusters which contain at least one 
Red [Super]Giant. For comparison, the solid black curve shows 
the empirical ZAMS of Schmidt-Kaler for stars. 



All-sky Compiled Catalogue of 2.5 million stars, 3rd version, CDS 
I/280B, ftp: //cdsarc.u-strasbg. f r/pub/cats/I/280B/ 



iKharchenko et al.l2005allbl) . The integrated magnitudes were ob- 
tained in the optical BV, and the NIR JHK S passbands by sum- 
ming up the individual luminosities of the most probable cluster 
members (see lPaperll for more details). In order to convert ap- 
parent magnitudes to absolute integrated ma gnitudes, we used 
the dist ances and extinction values given in Kharch enko et al.l 
( 2005a. b). The colours were computed as differences of the mag- 
nitudes in the corresponding passbands. 

From a comparison of our results with independent esti- 
mates of integrated magnitudes and colours available i n the liter- 
ature for Galactic open clusters (Batti nelli et al.|[l994t lLata et al.l 
2002), we found the rms-errors of the differences to be ±0.6 
mag and ±0.2 mag for magnitudes and colours, respectively 
(Paper I). These values may be regarded as the upper limits of 
the uncertainties of our data. The main reason for the differ- 
ences arises from the different membership criteria applied: in 
case of solely photometric selection of cluster members around 
the cluster centre (as usually happens in previous determina- 
tions), integrated luminosities can be overestimated due to ad- 
mixture of bright foreground stars, or underestimated due to 
missing of relatively bright members outside the central cluster 
region. However, in the following we assume our data to be accu- 
rate within ±0.6 mag for magnitudes and ±0.2 mag for colours 
though, one can expect that the realistic accuracy is a factor of 
two better. 

In Figs.Q]and[2]we show the combined colour-magnitude di- 
agrams (CMDs) for the most probable members of all our clus- 
ters, and the integrated CMDs of the clusters themselves. The 
figures explain our cluster classification which will be used for 
convenience throughout the paper. According to the evolution- 
ary status of their brightest stars, we consider three groups of 
clusters: the Main Sequence (MS) Clusters, the Red Giant (RG) 
Clusters and Post-MS Clusters. The Main Sequence Clusters 
contain Main Sequence stars only, the Red Giant group consists 
of clusters with at least one red giant or supergiant as a member 
in their upper CMD, and in the Post-MS clusters one observes 
that the brightest star just left the Main Sequence, but did not 
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yet reach the Red Giant/Supergiant domain. The classification 
of members into MS, Post-MS or giant stars is based on the data 
of the Padova isochrones. As separation lines we use the TAMS, 
and the line connecting the points of maximum luminosity at the 
shell hydrogen-burning stage. Since stars cross the Hertzsprung 
gap fast, the exact location of this limit is not so important. 

The other data of relevance for this paper are the cluster 
masses. They were estimated from the tidal radii of clusters 
which were drawn either from a direct fitting of King's pro- 
files to the observed density distribution of cluster members 
dPiskunov et alj 12007b or from calibrating the apparent semi- 
major axes o f the projected memb er distribution into the tidal 
radii system (iPiskunov et al.ll2008bl) . It must be emphasised that 
in our approach cluster masses and cluster luminosities are esti- 
mated independent of each other. This allows an unbiased con- 
struction of the M/L ratio. This circumstance overcomes the 
general disadvantage of the uncertainty of tidal masses of open 
clusters and makes their use in this study feasible. 

3. The construction of the SSP models 

A number of SSP-models has been constructed during the last 
decade. They are widely used fo r modelling both star clus- 
ters and galactic popula t ions (e.g. iMouhcine & Lanconl [2002; 
iBruzual & Charlotl [2003t lMarastonll2005l) . In this context we 
mention a study bv I Anders et alJ (|2009) who included the effect 
of mass loss via evaporation of star s from clusters in s t andard 
SSP-models, and investigations by Popescu & Hanson (2010) 
who take into account the discreteness of the initial mass func- 
tion (IMF). However, no SSP-models are available incorporat- 
ing both effects. Since these issues seem to be of great relevance 
for model ling Gal actic open clusters, we refine the model pre- 
sented in iPaper III bv considering both effects. In order to make 
sure that the models we produce are compatible with present-day 
SSP models, we construct models both in the "standard" and 
"extended" modes. The standard mode reproduces the features 
of the adopted SSP models, whereas the extended mode takes 
into account both the effect of a discrete IMF and the effect of 
mass loss. In this section we describe the model ingredients and 
compare the results of our modelling with published standard 
SSP-models. 

3.1. Star formation 

The basic model parameter is the cluster mass M c . We assume 
that at the moment of cluster formation (as a classical open clus- 
ter) the placental gas has already been removed and cluster con- 
sists of stars only. For simplicity, we also assume a simultaneous 
formation of stars of different masses. The IMF of the cluster 
stars formed with mass to 

/(to) = dN/dm 

is defined within the lower and upper limits (mi and m u ) of the 
allowed masses of cluster stars. Here dN is the number of stars 
having masses between m,m + dm. We adopt a power-law form 
for the IMF, 

f(m) = km- a (1) 

with parameters a (IMF slope), and k (normalisation factor) 
where k using proper normalisation can be computed via cluster 
mass M c , slope a, and the limits to; and m u 

(2 - a) M c 
ml- a - m 2 ~ a ' 



There are several reasons of choosing a simple power-law IMF 
valid over the whole range of stellar masses. First of all, the in- 
tegrated fluxes of open clusters do not depend strongly on the 
luminosity nor on the number of low-mass stars (i.e. IMF pa- 
rameters at lower masses). Secondly, the mass spectrum at the 
low-mass end is changing with time due to dynamical evolu- 
tion, so the initial distribution of stars is no longer valid for older 
clusters. Finally, the low-mass IMF is poorly known, and dif- 
ferent poss ibilities (including t he Salpeter IMF) are still under 
discussion (Malkov et al . 1997). Also, recent observations indi- 
cate that the cluster mass spectru m follows a power-law shape 
with an exponent of about —1.8 (Bou dreault et alJfeOlOh down 
to small masses of to ~ 0.1 to©. Therefore, we decide to restrict 
the modelling to a power-law IMF confined with some "effec- 
tive" low mass limit mi. Throughout the paper, the mass spec- 
trum is represented by the Salpeter IMF (a = 2.35). 

The largest possible mass of cluster stars m u is defined by 
the adopted grid of evolutionary tracks/isochrones. Only at the 
very early stages of cluster evolution does m u have an impact 
on the integrated colours. In contrast, the lower mass limit mi 
governs the population at the massive end of the IMF provided 
that the IMF slope and the initial cluster mass M c are fixed. In 
the following, we consider M c and to; to be free parameters of 
the adopted model. 

3.2. Stellar evolution 

The photometric properties of model clusters are defined by the 
implemented grid of iso chrones. Here we u s ed the grid provided 
by the Padova group dBertelli et al.l | l994t iGirardiet all 120001: 
Marigo & Girardi 2007; Marigo et al. 2008) via the online server 
CMEfl According to our experience, the Padova isochrones fit 
the brighter parts of the open cluster CMDs reasonably well over 
a wide range of cluster ages. These parts are responsible for the 
integrated luminosities and colours. We used a solar metallicity 
grid (Z — 0.019), and retrieved the passbands U, B, V, R, I and 
J, H, K s with an age range logt = 6.0 . . . 10.2 and a step of 
0.02 in log t. 

3.3. Discretization of the model 

Adopting the continuity of the IMF in standard SSP models, one 
assumes that every (infinitely) small mass interval of the given 
isochrone emits light according to its T e ff and log 5. The total 
flux that a "theoretical" cluster emits is a result of the integration 
over the isochrone, weighted by the given IMF. It is immediately 
obvious that in this case luminosity must scale with mass, once 
the age (or the isochrone) is fixed. In the discrete approach, how- 
ever, only those loci on the isochrone contribute where stars of a 
respective mass are actually present. 

With the distribution function of stars given by the IMF, we 
construct an initial sample of cluster members for a given set of 
model parameters. The total number of members iV c is defined 
by the integral 

N c — J /(to) dm, 

mi 

and the mass of the fc-th member to^ is assigned at random by 
the distribution given by eq. (Q]i. At a given age of a cluster, each 
cluster member contributes to the cluster luminosity in accor- 
dance with its evolutionary state. 

2 http://stev. oapd . inaf . it /cgi-bin/cmd 
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Figure 3. Mass loss from star clusters according to high- 
resolution N-body models by IBerczik & Petrovl d2008l) . The left 
panel shows the relative loss of the total mass of clusters or- 
biting at the Solar Galactocentric distance {Rq = 8.5 kpc). 
The solid (green), dotted (blue), and dashed (red) curves are for 
clusters of initial masses M c = 10 3 , 5 • 10 3 , 10 4 TOq, respec- 
tively. The right panel demonstrates selective mass loss occur- 
ring among stars of different masses to/to© in a cluster with 
M c = 10 3 toq. The curves show the evolution of the MF-to- 
IMF ratio F(m,t)/ /(m) with time, starting from t = Myr, 
with a step of 100 Myr and las t point t = 1249 Myr. The dots 
show the data of lAnders et al.l Q2009) scaled to a cluster with 
mass M c — 10 3 to and the age of 1249 Myr. The red portions 
of the curves correspond to stars alive and grey ones to stellar 
remnants. 



For a given M c , open clusters represent random realisations 
of the IMF which have statistical fluctuations even for a smooth 
probability density distribution. At the high-mass end of the 
mass range, these fluctuations are particularly large because of 
low number statistics, and produce gaps in the mass spectrum. 
The fluctuations depend not only on the model mass, but also 
on the lower mass limit which, at a given M c , controls both the 
upper mass limit and number of massive stars. For example, for 
the Salpeter IMF, M c = 1000 m and to/ = 0.01, 0.1, and 
1 TOq, the upper mass limit m u turns out to be about 60, 90, and 
100 toq, on average. Note that m u is confined to 100 to q due 
to the adopted evolutionary grid which does not include higher 
mass stars. 

3.4. Mass loss from the cluster 

Usually the standard SSP-models neglect mass loss from star 
clusters. Indeed, the bulk of the cluster luminosity comes from 
the brightest stars. Their mass distribution is relatively weakly 
affected by dynamical processes if compared to low-mass stars. 
On the other hand, masses of extragalactic clusters are usu- 
ally determined from the mass-luminosity relation (M/L ra- 
tio) which is also based on the standard SSP models. A recent 
attempt to inc orporate dynamical processes into SSP models 
was made by lAnders et al.l d2009|) who include d results of N- 
body calculation of Bau mgardt & Makinol (12003) into a standard 
GALEV model (i.e., with continuous IMF). 

In the present paper, we use th e results of N-body simu- 
lations bv IBerczik & Pe trov (2008) which were carried out to 
reproduce the local population of the Galactic open cluster^. 

3 The full model snapshots are available from 
ftp://ftp.ari.uni-heidelberg.de/pub/staff/berczik/phi-GRAPE-cluster 



Table 1. Parameters of adopted N-body models 



M c (0), M Q 


N c (0) 


-R c (0),pc 


t dis , Myr 


1000 


4040 


3 


1511 


5000 


20202 


7 


4661 


10000 


40404 


10 


7572 



Below we give a short description of the model parameters. For 
more details see lKharchenko et al.ld2009l) . The models represent 
ensembles of stars distributed at birth along with the Salpeter 
IMF (a = 2.35, mi = 0.08 Mq, m u = 8M Q ) and mov- 
ing in the external tidal field of the Milky Way at the Solar 
Galactocentric circular orbit (Rq = 8.5 kpc). 

Such an approach means that both dynamical- and SSP- 
models are coordinated with respect to the properties of the stel- 
lar population. The adopted set of N-body model parameters is 
listed in Table Q] In the first three columns of the table we show 
initial mass, number of stars comprising a model, radius of the 
model clusters, and the last column indicates their estimated dis- 
solution time. Further, we assume that th e phase of the initia l 
mass loss ("infant weight loss" according to Wei dner et al.l2007h 
is relatively short and is finished with the death of the most mas- 
sive cluster members (m > 8 Mq). We trace the evolution of the 
N-body model within a time interval of 1 .2 Gyr. Consequently, 
in this study we consider the intermediate phases of cluster dy- 
namical evolution occurring between the infant weight loss and 
the final decay phases. However, this embraces more than 90% 
of the observed local cluster population. 

In the context of the present study, the most interesting fea- 
ture of the dynamical models is the steady mass loss appropriate 
to clusters due to their dynamical evolution and its consequence 
for the stellar content of clusters. These functions are shown in 
Fig. [3] In the left panel we compare the relative change of the 
total mass of the model during the first 1.2 Gyr of the evolution. 
Both effects of stellar and dynamical evolution are taken into 
account. It is assumed that the gas lost by the stars is removed 
from the system, whereas the stellar remnants remain within the 
cluster and continue their dynamical evolution as normal stars of 
the respective mass. The dynamical mass loss has a regular na- 
ture and occurs due to two-body interactions within the cluster. 
No such events as molecular cloud encounters are taken into ac- 
count, so our results provide lower limits of possible dynamical 
effects. According to Fig. [3] the rate of the relative mass loss is 
higher in massive clusters at the beginning of their evolution, but 
after about 600 Myr it slows down, whereas the low-mass clus- 
ters continue to loose mass with approximately constant rate. 
At log t « 9.2 the low-mass clusters have almost completely de- 
cayed. For more massive clusters this process takes considerably 
longer time (see TableQ}. 

The right panel of Fig. [3] shows how effectively stars of dif- 
ferent masses evaporate. As a measure of the dissipation process, 
we use 

q{m,t)=F{m,t)/f{m), (2) 

where /(to) is the initial mass function, and F(m, t) is the mass 
function (MF) of cluster stars at an instant time t. The MF is 
a result of both stellar and dynamical evolution of the cluster 
population. The former process results in shortening the mas- 
sive end of F(m, t) and in decreasing q(m, t), whereas the dy- 
namical evolution is more effective in depleting stars at the low- 
mass end of F(m, t). The input parameters of our SSP models 
include the possibility of stellar masses outside the mass range 
of 0.08 Mq. . . 8 Mq adopted for N-body simulations. In these 
cases, the constructed relations are linearly extrapolated. As the 
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evaporation is more efficient than the expected value from linear 
extrapolation, our estimates of q(m, t) are too conservative for 
stars with masses m < 0.08 m©. In contrast, this extrapolation 
is of minor importance for high-mass stars, since the lifetime of 
stars with m > 8 m© is much shorter than the time of a signifi- 
cant decline of q(m, t) from unity. 

For comparison, we plot the results of Ande rs et al.l J2009) 
in Fig. [3] as filled dots. In order to place their dat a into F ig. [3j 
we co mpute the MF by using formula (6) from lAnders et al.l 
(2009). A normalisation (i.e. the amount of cluster mass loss) 
is achieved by shifting the resulting function q along the vertical 
axis to provide the best fit at t = 1.2 Gyr. In general, the results 
from both appro aches are in reasonab ly good agreement, except 
for a jump in the lAnders et al.l d2009t) data which occurs due to a 
rather rough separation of the data into two mass intervals. Thus, 
our approach allows a more accurate approximation of the mass 
loss process along stellar masses. 

As a "by-product" of the N-body simulations we get data 
allowing to follow the change of the MF slope caused by the 
dynamical evolution of a cluster alone. Assuming the Salpeter 
IMF (a = 2.35), we obtain a — 2.0 i 0.1 in the mass do- 
main (0.08, 1.0) TOq at logf = 9.5. This can be compared 
with the recent finding of a slope of a = 1.8 ± 0.1 derived by 
Boudreau lt et al.l (12010) for the Praesepe mass function based on 
star counts in the same mass domain. The authors explain the flat 
MF by the possibility that the Praesepe IMF is different from the 
Salpeter one. The N-body simulations show, however, that about 
the same flattening of the Salpeter IMF can be achieved if dy- 
namical evolution of clusters is taken into account. 



3.5. Integrated flux calculation 

For the case of a continuous IMF and simultaneous star forma- 
tion, the integrated luminosity of a cluster in a passband i at a 
moment t is computed as 

m(t) 

L>i(t) = / F(m,t) li(m,t) dm , 



where F(m, t) is the PDMF, li(m t t) is the mass-luminosity re- 
lation for stars in the ith passband valid at the moment t, and 
m(t) is the mass of stars just finishing their life at the moment 
t. Both functions are taken from the theoretical isochrones de- 
scribed above. 

In the case of a discrete IMF, we simply reproduce the cal- 
culation of integrated fluxes of real clusters by adding the lumi- 
nosities of non-evaporated individual stars 



Li(t) 



N c (t) 

£ 

fc=i 



k(m kl t) , 



the stars are labelled with their initial masses m k , and N c (t) is 
the number of cluster members with lifetimes exceeding t. 
The integrated magnitudes I (Mi) are computed as 

I{Mi) = -2.5 log U-Mf , 

with Mf to be the absolute magnitude of the Sun in the pass- 
band i. In the case of two passbands (i and j), the integrated 
colour I(i — j) is determined as 

= -2.5 log Li/Lj. 




Figure 4. Evolutionary tracks developed in this paper (solid 
red curves) and models computed with help of the online 
SSP-servers SB99 (dotted green curve) and GALEV (dash- 
dotted blue curve). The left panels show optical magnitudes and 
colours, the right panels are for NIR passbands. The adopted 
model parameters are given in the text. 



3.6. Comparison with literature models 

In Fig. 3] we compare our discrete model with the data com- 
puted with the help of online servers provided by the GALEVQ, 
and Starburst993(SB99) sites which de facto can be regarded as 
present-day standards. Our model was calculated for a continu- 
ous IMF for a cluster with M c — 10 3 m©, and stellar mass limits 
mi = 0.1 m© and m u = 93 m©. The mass spectrum is repre- 
sented by the Salpeter IMF (a = 2.35). The on-line continuous 
models were computed for M c = lO 6 ™©, the Salpeter IMF, 
and stellar mass limits of mi = 0.1m© and m u = 100 m©, 
which are compatible with our values. We consider the photo- 
metric passbands B, V (Johnson) and JHK S (2MASS). Due 
to the difference in cluster masses between our and the on-line 
models, we simply moved the on-line tracks along the magni- 
tude axis to reach the best agreement between the curves. Since 
the integrated colours do not depend on cluster masses in the 
case of continuous IMF, the tracks in the colour-age diagram are 
plotted without any shift. 

According to Fig. |4] (left panels), the models coincide well 
in the optical. A small disagreement between our and the SB99 
tracks at early ages can be explained by slightly different up- 
per mass limits adopted. In the NIR the agreement between the 
three models is worse. Especially, the SB99 models differ promi- 
nently from both GALEV and our results. This can occur due 
to a different grid of the underlying isochrones used by SB99. 
Smaller differences between our models and GALEV probably 
come from discrepant transformations into the 2MASS system 
applied in GALEV and our models. A general impression from 
Fig.|4]is that our models are located in between the data of SB99 
and GALEV. Therefore, we conclude that our models describe 
the same system as the other present-day SSPs. 



http : / /www . galev . org/ 



http : / /www . st sci . edu/ science/ starburst 99/ 
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Figure 5. The evolution of cluster luminosity in the optical (upper row) and NIR (lower row). The curves in the left and middle panels 
present SSP models with continuous IMF, whereas in the right panel we compare continuous and discrete options. Left panels: 
impact of dynamical mass loss on integrated luminosity of clusters of different masses and stellar low-mass limit mi — 0.1 m . 
The four curves from the top to the bottom are models for clusters with masses M c — 10 m , M c = 5 ■ 10 3 m , M c = 10 3 m , 
and M c — 5 • 1O 2 to , respectively. The dotted curves indicate models without mass loss. Middle panels: influence of different 
lower mass limits m; of cluster stars on integrated luminosity of clusters with mass M c = 1O 3 to . The upper, middle, and 
bottom curves correspond to mi = 1, 0.1, 0.01 m , respectively. Right panels: integrated magnitudes of clusters modelled with 
continuous (curves) and discrete (connected symbols) IMFs for stellar low-mass limit mi = 0.1 m . The upper curve and crosses 
are for cluster mass M c — 10 4 m , and the bottom curve and diamonds are for M c = 10 3 m . The mass loss due to dynamical 
evolution is incorporated in both models. For comparison, the continuous model without mass loss is shown by the dotted curves. 



4. Luminosity evolution of star clusters 

Integrated luminosities (magnitudes) are the basic data provid- 
ing estimates of masses of extragalactic clusters. They are, also, 
used to describe the cluster population as a whole, e.g. the his- 
tory of cluster formation, and its spatial variations (via the lumi- 
nosity or mass functions). It is generally accepted that the current 
SSP-models reproduce well the integrated magnitudes of extra- 
galactic clusters. However, there is a number of effects which 
must be taken into account to avoid serious aberrations by the 
interpretation of observations. 

4. 1. Effects in the magnitude vs. age diagram 

In Fig. [5] (left panels) we check the impact of the mass loss ef- 
fect on integrated magnitudes of clusters of different masses. 
We consider models based on the continuous IMF. Within the 
range of cluster masses from 5 • 10 2 m to 10 4 m - typical 
for Galactic open clusters - the spread of cluster luminosities 
reaches up to three magnitudes both in the optical, and in the 
NIR. The mass loss due to dynamical evolution leads to some 
decrease of the integrated luminosity provoked by a relatively 
minor change of the upper mass spectrum (if compared to the 
low-mass portion of the IMF). The more massive a cluster is, the 



earlier the mass loss effect becomes significant. However, the ef- 
fect is small, from AM w 0.2 for M c < 10 3 M to AM w 0.5 
for M c = 1O 4 M . This can be compared with the decrease 
of integrated magnitude due to stellar evolution which reaches 
about 5 to 6 mag in the optical and about 3 to 5 mag in the NIR. 

The impact of the lower mass limit mi in the mass spec- 
trum of a cluster is shown in the middle panels of Fig. [5] It may 
seem to be strange that the population of the low-mass end of the 
mass function can affect the cluster integrated luminosity which 
is mainly defined by the most massive stars. For the given IMF 
and cluster mass, however, the lower mass limit governs the po- 
sition of the "mass centre" of the system and hence, the pop- 
ulation of the up per part of the mass spectrum, too. For other 
IMF shapes (sav lMiller & Scalolll979h . the low mass limit can 
be translated into the effective mass limit where the IMF is at its 
maximum. According to Fig. [3] the effect produces variations of 
AMy w 2 mag when mi varies between 0.01 m and 1 m , 
and it has the same effect both in the optical and in the NIR. 

In the right panels of Fig. [5] we compare the influence of the 
continuous and discrete IMFs on the integrated magnitudes of 
open clusters of different masses, and the same m;. In the op- 
tical, one can observe a reasonable agreement over the whole 
age range considered. As expected, the discrete model produces 
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Figure 6. Integrated magnitudes of open clus- 
ters in two different passbands vs. cluster age. 
The observations of the Galactic open clus- 
ters are shown by dots, blue for MS-clusters, 
green for Post-MS clusters, and red for RG- 
clusters (see Sect. [2] for explanation). Bars in 
the upper-right corner of each panel indicate 
typical errors of the integrated magnitudes. The 
curves are evolutionary tracks of a continu- 
ous model with effects of mass loss incorpo- 
rated for model parameters M c — lO 3 m0, 
and mi = 0.01 m© (solid curve), and M c = 
5 • 10 2 ttiq, mi =0.1 m.0 (dash-dotted curve). 
The latter option (M c = 5 • lO 2 m0, m; = 
0.1 mo), but with discrete IMF, is shown by 
the connected crosses. 



a spread around the continuous track. At large M c , the discrete 
option converges to the continuous curve, and at M c — 10 6 the 
discrete and continuous models completely coincide. However, 
at smaller M c one observes a systematic difference between in- 
tegrated magnitudes modelled with the continuous or the dis- 
crete IMF. At short time scales, the track of the discrete model is 
represented by recurrent events consisting of a slow increase of 
cluster luminosity due to Main Sequence evolution of the bright- 
est star and of a relatively fast outburst caused by the formation 
of a bright short-living supergiant. After the supergiant's decline, 
the process repeats caused by the follow-up brightest star on the 
Main Sequence. With increasing cluster age, the outbursts be- 
come wider and less intensive due to the lower masses of the 
stars involved. 

This scenario is valid for all passbands, but the effects be- 
come stronger with increasing wavelengths. The systematics are 
especially prominent in K s (see the lower right panel of Fig. [5]). 
Depending on mass and age of clusters, the shift between the 
models with continuous and discrete IMF reaches from 1 to 
3 mag. For younger clusters (\ogt < 7.5), one is not able to 
discriminate between the discrete models for M c = 10 3 and 
10 4 rri0 in Fig. [5] We attribute this bias to the considerable in- 
crease of the luminosity of red supergiants in the NIR with re- 
spect to coeval MS-stars. 

4.2. Comparison with observations 

In Fig. [6] we compare the observed magnitude-age relations ob- 
tained in different passbands for 650 Galactic open clusters with 
the corresponding models. The typical accuracy of the derived 
magnitudes of about 0.6 mag is plotted in the upper right corner 
of each panel. The observed spread of the data points consid- 
erably exceeds this value. The obvious reason is a difference in 
the initial mass of the observed clusters. However, as we will 
discuss below, this might not be the sole effect responsible for 
the scattering. The blue dots indicate the MS-clusters outlin- 
ing the lower border of the distribution, whereas the red dots 
(RG-clusters) tend to occupy the upper stripe, and especially the 
upper-right part of the distribution. This distribution reflects evo- 
lutionary changes in the cluster population: the older a cluster is 
the higher is the probability to observe a red (super)giant among 
its brightest members. At redder passbands, the distribution be- 
comes more dishevelled though, the segregation between MS- 



and RG-clusters gets more evident. This points to the relatively 
higher luminosity of red giants in the NIR. 

The observed diagram can be fitted by a set of models with 
different cluster parameters. In general, the upper bound of the 
distribution is compatible with cluster masses of the order of 
lO 4 m0, whereas the lower bound can be described by mod- 
els with a few hundreds of solar masses. In the middle, the 
cluster models for about 1000 toq can fit the distribution. In 
Fig. [6] we include the models with continuous IMF which fit 
the middle part of the observed magnitude-age relations. The 
model parameters chosen are (M e , mi) — (5 ■ 10 2 , 0.1) m Q and 
(10 3 , 0.01) 77i0. An additional scattering of the modelled points 
can be produced for clusters with M c — 5 • 10 2 m© by varying 
mi within the extreme range from 0.01 to 1 m . However, the 
large spread of observed integrated magnitudes can be, in a more 
natural way, described by the models with a discrete IMF. 

This can be especially well seen at young ages (log t < 8) 
where open clusters of the same age and the same mass (e.g., 
M c = 5 • lO 2 m0) can have integrated luminosities differing 
by about 2 mag in the optical and about 5 mag in the NIR. The 
MS-models with the discrete IMF are systematically located be- 
low the corresponding continuous models, but they coincide well 
with observations of the MS-clusters. During RG-events, the dis- 
crete models achieve considerably larger luminosities than the 
values predicted by the standard SSP models but, again, compat- 
ible with observed integrated magnitudes of RG-clusters. This 
has important consequences for the mass determination of open 
clusters when based on standard SSP models: even perfectly 
determined integrated magnitudes can provide incorrect cluster 
masses. Since RG-events are relatively rare and short at earlier 
stages of a cluster's life (about 20 events during the first 100 
Myr for a cluster with M C) mi = 10 3 , 0.08 m©), the mass of 
this cluster will be, with a higher probability, underestimated by 
standard SSP models due to the magnitude bias. In the case of 
RG-clusters, the masses would be considerably overestimated. 
Again, this conclusion is valid for typical Galactic clusters in 
the solar neighbourhood. At high cluster masses the discrete- 
ness effect gradually vanishes. However, we must keep in mind 
that the integrated magnitude of a cluster of moderate mass (e.g. 
« 10 3 m Q ) increases significantly during RG-events. So, it can 
be more easily detected in the NIR in other galaxies and assigned 
an overestimated mass from the standard SSP model. 
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Figure 7. The evolution of cluster colour in the optical regime 
and the effects of M c (upper row) and to; (bottom row). The 
curve in all panels shows the continuous-IMF model, the con- 
nected symbols show the discrete-IMF case. In the upper row we 
compare M c = 5 ■ 10 2 to© (left) and 10 4 to© (right) cases for 



m; = 0.1 to©. The lower row is constructed for M c 
and to; = 0.01 to© (left) and 1 m© (right). 
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5. Colour evolution of star clusters 

In iPaper Ij, we have already reported on the consequences of 
using standard SSP models for estimates of cluster ages from 
observed integrated colours when the effect of a discrete IMF 
is neglected. In this paper we extend this analysis and consider 
all the effects studied in the previous sections (mass loss from 
a cluster, the IMF-discreteness, and the effect of the low-mass 
limit) and investigate their impact on cluster colours. 

5. 1. Effects in the colour-age diagram 

In Figs. [7] and [8] we compare the impact of cluster mass and of 
the lower limit of stellar masses onto the evolution of integrated 
colours of a cluster in the continuous and discrete approxima- 
tions. Again, we consider both optical (Fig. [7J, and NIR (Fig. [8) 
colours. As expected, the continuous model does neither depend 
on M c nor on to; . The only parameter influencing the standard 
model is the IMF slope a, although this influence is not strong. 
To substantially change the track, one would need to increase 
the exponent a to 4.35. Also, a flattening of the IMF with re- 
spect to the Salpeter value of the slope does not have a strong 
effect on the colour-age relation. On the other hand, the discrete 
model demonstrates a strong dependence on the value of M c or 
to; and shows a considerable divergence from the model with 
the continuous IMF. 

First of all, we observe a systematic bias produced due to 
the discreteness effect. At small masses M c red giants are rare, 
and the clusters spend a relatively long time as Main Sequence 
clusters free of red giants (similar to the present-day Pleiades). 
In contrast, the continuous models cannot resemble MS clusters 
(and the Pleiades) since they always contain the RG-portions of 
the isochrones. This produces a bias which is especially well 
seen in Fig.[7j at M c — 5 • 10 2 to© (in the upper-left panel) where 
we practically observe two different sequences corresponding to 
the continuous and discrete options. From time-to-time the dis- 




Figure 8. The evolution of cluster colour in the NIR and the ef- 
fects of M c (upper row) and to; (bottom row). The model pa- 
rameters and designations are the same as in Fig [7] 



crete track is interrupted by a "red-giant event", when a clus- 
ter becomes extremely red. Simultaneously, as we did see in 
the previous section, the cluster peaks in its luminosity, so we 
call such an event a "red flash". At log t > 9, when the clus- 
ter red giants become more long-living and the travelling dis- 
tance from the Main Sequence to the Rad Giant branch becomes 
shorter (i.e., the Hertzsprung gap becomes narrower), the se- 
quences converge, and the colour fluctuations diminish. 

With increasing cluster mass the existing holes in the mass 
spectrum are filled with stars, the IMF approaches a continu- 
ous shape, the cluster track approximat es the co ntinuous option, 
and the colour bias is vanishing (see Paper II, for a more de- 
tailed discussion of this effect). However, the upper-right panel 
of Fig. |7]shows that at M c = 10 4 m© - the upper limit of cluster 
masses observed in the Solar vicinity - a sufficient colour bias 
of the order of A(B — V) = 0.1 mag is still seen and accompa- 
nied by appreciable scatter in colours. Actually, the fluctuations 
of the colour track reflect the degree of discreteness of the IMF. 

A comparison of the upper and lower rows of Fig. [7] also 
shows that both M c - and m;-effects have approximately the 
same strengths. The bottom-left plot indicates that for a high- 
mass cluster the decrease of to; partly compensates (with re- 
spect to the colour bias and fluctuations) the increase of M c . One 
should note, however, that the colour fluctuations in this case are 
somewhat smaller, and cease earlier than in the low-mass model 
with "normal" m; = 0.1 to©. 

The flashing activity depends on cluster age, but this depen- 
dence is different for different cluster masses. Massive flash- 
ing clusters are better represented at younger ages (logt < 
7.8), whereas low-mass clusters show flashes typically later (cf. 
Fig. 17). However, the low frequency of "red-giant events" in a 
low-mass cluster is compensated by the large number of low- 
mass clusters due to the increase of the cluster initial mass func- 
tion to lower M c . Therefore, at younger ages one observes both 
massive and low-mass flashing clusters. At log t > 7.8 the mas- 
sive clusters approach the continuous track, and the flashing pop- 
ulation is represented exclusively by low-mass objects. 

As for integrated magnitudes, the effects for integrated 
colours become stronger at longer wavelengths. This statement 
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Figure 9. Integrated colours of open clusters vs. cluster age. The 
left panel shows data for (B — V), right one for (J — K s ). 
The observations of the Galactic open clusters are indicated 
by dots, blue for MS-clusters, green for Post-MS clusters, and 
red for RG-clusters (see Sect. [2] for explanation). Bars in the 
upper-right corner of each panel mark typical errors of the inte- 
grated colours. The solid curves are evolutionary tracks of a SSP 
model with a continuous IMF for a cluster with M c = 10 3 m© 
and mi = 0.1 m©. The connected symbols (crosses) show the 
discrete-IMF option of this model. 

is supported by Fig. [8] the counterpart of Fig. [7] for (J — K s ) 
colour. In the NIR we observe strong fluctuations of the order of 
1 mag, and a bias of about a few tenths or more even in models 
of M c — 10 4 to©. The models with discrete IMF always include 
the effect of mass loss. In the case of a continuous IMF, this ef- 
fect turns out to be negligible, as the tracks with or without mass 
loss coincide, in practice. 

5.2. Comparison with observations 

In Fig.|9]we compare the observed colours of the Galactic clus- 
ters with the adequate models in the colour-age diagram. The 
observed clusters with ages log t < 8.5 are predominantly repre- 
sented by MS-clusters which show a relatively narrow distribu- 
tion in colour. A redward spread of the data points in the diagram 
is caused by a few clusters containing yellow or red supergiants, 
and it is considerably exceeding typical photometric errors. 

At log t > 8.5, one finds, in general, clusters with red giants. 
The scatter in the observed colours of these clusters is compa- 
rable with what is observed for younger ages (opposite to the 
predictions by the model with M c = 10 3 to©). Part of this scat- 
ter can be explained by the different evolutionary stage of red 
giants along the RG-branch in different low-mass clusters. The 
other part follows from the lower accuracy of their observations. 
The older clusters are, on average, absolutely fainter than the 
younger ones. At distances typical for our sample (< lkpc), they 
have fainter apparent magnitudes, sometimes close to the limit of 
the accuracy of our data. As a consequence, the accuracy of the 
cluster parameters is typically lower for older clusters than for 
younger ones. Especially, the lower accuracy of the reddening 
determination contributes to the colour spread observed for older 
clusters. The depletion of the cluster sequence after logt ~ 8.5 
can be explained by either cluster dissolution (the typical life- 
time of local cluster s is only around logt w 8.5, according to 
iPiskunov et al.ll2006h or incompleteness effects, which are more 
significant for older clusters. 

Again, the model with a continuous IMF is shown in the 
plots as the solid curve, whereas the discrete option is repre- 
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Figure 10. The observed mass-luminosity ratio vs. age of the 
Galactic open clusters in the optical (V) and NIR (K s ) photo- 
metric passbands. The dots are data for individual clusters, the 
bars indicate the individual errors in the data. The open dots 
mark extended clusters excluded from the analysis (see text). 
The large black dots show the running median of the M/L ratio 
computed with a (log i)-bin of 0.5 and a step of 0. 25. The straight 
line in the left panel is the empirical relation from Bat tinelli et al.l 
(11994 . 

sented by connected symbols. There is no way for the contin- 
uous model to fit the observed colours of either MS-clusters 
or red "outliers", undergoing the red flash. In contrast, the dis- 
crete model is able to describe the observed distribution of both 
groups. The colour scatter of MS-clusters probably reflects the 
true error of cluster colour determination, as well as the secular 
colour evolution of the Main-Sequence stars. 

In Fig.|9]we show the discrete model with M c = 10 3 to© and 
in 1 = 0.1 to©. It produces a mass spectrum consisting of 2850 
stars. In the mass range 5 m a > m > 0.3 to©, however, the 
number of stars is only 625. This is consiste nt with 780 Pleiade s 
members observed in this mass range bv lBelikov et' al. (1998). 
Thus, the model appears to be comparable to the Pleiades. 
According to the IMF normalisation, the number of model stars 
with to > 5 to© turns out to be 17. So, before the model cluster 
reaches the age of the Pleiades (\ogt « 8.1), it has produced 
17 red supergiants. The comparable number of corresponding 
"RG-events" is clearly seen in Fig. [9] Outside these short events, 
the cluster is seen as a main-sequence cluster during the first 
125 Myr and resides in the MS-domain of the diagram. If we 
increase the lower mass limit mi, the number of massive stars 
increases and the fraction of RG-events slowly increases, espe- 
cially at older ages. But even at To; = 1 to©, the MS-clusters 
are present in the track up to log t w 8.2. On the other hand, de- 
creasing to; diminishes the number of RG-events. In principle, 
from the frequency of red outliers in the "colour versus log t " di- 
agram, one can estimate indirectly the parameter to; (see Sect. [7] 
more detailed discussion of this issue). 

For the other colours the behaviour is qualitatively similar, 
although the redder the passband, the stronger is the disagree- 
ment between continuous and discrete SSPs. As Fig.|9]shows in 
the (J — K s ) panel, the colour offset between the prediction 
of the continuous-IMF model and observations reaches about 
1 mag at log t < 8. At older ages it tends to be smaller, but still 
has a significant value of about 0.5 mag. This occurs because 
the red portions of the isochrone radiate more of their light at 
longer wavelengths than the blue portions, and the inherent bias 
increases. Similar to the integrated magnitudes, the higher the 
cluster mass the smaller is the bias. 
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The observed dispersion in the colour distribution of clus- 
ters at \ogt < 8.5 in Fig. [9] can be understood if one as- 
sumes the presence of a scatter of initial masses among the ob- 
served clusters. The models for low- and medium-mass clusters 
(M c < 10 3 m©) provide the MS-clusters and red-outliers. The 
models for more massive clusters (with M c about and more than 
10 4 ra 9 ) provide moderately red clusters in the vicinity of the 
continuous track. At older ages the tracks with continuous and 
discrete IMF converge. According to Fig.|9j we conclude that the 
Galactic open clusters build a sequence of MS-clusters and do 
not show a concentration towards the model with a continuous 
IMF. This agrees qualitati vely with a power-law mass distribu- 
tion of the clusters (see e.g. lPiskunov et all20 08a) which implies 
the dominance of low-mass clusters in the Solar neighbourhood. 



6. Mass/luminosity ratio - age relation 

The mass-luminosity ratio M/L plays an extremely important 
role in studies of extragalactic clusters. Unfortunately, it is es- 
tablished empirically for massive star clusters like globulars and 
almost completely unknown for less massive clusters. The ma- 
jor reason has been the low number of reliable masses of open 
clusters. That is why the mass-luminosity ratio is normally taken 
either from counted masses or from model predictions. The mod- 
els themselves need to be constrained by the observations. In this 
section, we construct the observed M/L-age relations based on 
tidal masses of open clusters and use them for constraining SSP- 
models, where we take into account several effects usually ne- 
glected as model ingredients. 

6.1. Empirical relation 

The sample used in this section consists of clusters with 
sufficiently accurate masses (error 8m s less than 80%) from 
iPiskunov etaD (l2007l) : IPiskunov etail (l2008bl) . Further, we ex- 
clude nine objects, frequently classified as compact associations 
(Cep OB3, Cyg OB2, NGC 869,NGC 884, Nor OB5, Sco OB4, 
Sco OB5, Sgr OB7, Vel OB2) from consideration. The distribu- 
tion of the remaining clusters in the M/L-age diagram is shown 
in Fig. [10] A few of these clusters seem to have exaggerated val- 
ues of the Mj L ratio and form in Fig. [TTJa separate group of 1 1 
objects with Mj Ly > 2, typical of globular clusters. This group 
consists of nearby clusters most of which are poor and irregular 
(Cr 65, Platais 6, Platais 9, and Turner 5), but also includes the 
a Per cluster, which can be regarded as an association, and two 
rich clusters Blanco 1, and Stock 2. We can not exclude that 
we overestimate their tidal radii, and hence their tidal masses. 
Therefore we mark this group with open circles in Figs. [TOl and 
[T2l and do not consider them in the following discussion. The 
remaining sample includes 493 clusters altogether. 

The straight line in the M/Ly-age diagram shows the em- 
pirical relation con structed from a sample of 138 clusters by 
Batti nelli et al.l (Il994l) . They used stellar counts for evaluating 
the cluster mass. Therefore, the derived masses do not take into 
account unseen stars and can be regarded as lower cluster mass 
estimates only. On the other hand, the cluster luminosity is de- 
fined by the brightest stars. It is clear that such an approach 
provides only a lower limit for the Af/L-ratio. This is well- 
illustrated by Fig.fTOl 

The thick curves connecting the large black dots in Fig. [10] 
represent the running median of the Mj L ratio with a (log t)-b\n 
of 0.5 and a step of 0.25, and it is constructed on the basis of the 
sample of 493 Galactic open clusters described above. FigurefTOl 



shows that in the considered range of ages the M/Ly-age rela- 
tion is a rather shapeless cloud of data points. However, a certain 
trend can be observed when we take the medians in age bins. We 
find a formal average value of (M/Ly) = 0.26 ± 0.01 with a 
standard deviation of ±0.30. Adopting the error of an integrated 
magnitude to be ±0.6 (Sect. |2]l and taking into account that the 
error in M c is better than ±80%, we derive the relative error 
of a single measurement of the ratio to be better than ±0.97. 
Assuming further that the spread of points in Fig. [TTJaround the 
"true" value of 0.26 is caused solely by measurement errors, we 
expect a standard deviation of 0.25. This value is slightly lower 
than the observed scatter, and the additional dispersion produced 
by unknown effects should be less than 0.2. 

The relation in the NIR is in many respects similar to the op- 
tical one. It also shows a considerable dispersion but there is a 
tendency for M/Lk, to decrease, as the cluster age increases. 
The running median of the ratio supports this impression. The 
formal average value of the M/Lk s ratio is 0.66, and the stan- 
dard deviation is 1.07. The expected dispersion due to measure- 
ment errors turns out to be 0.64. This leaves much room for un- 
known effects (of the order of 0.86 mag). 

In the sections below, we consider different SSP-models, try 
to fit them to the observed distribution of the M/L ratio and 
discuss possible reasons for its relatively large scatter. 

6.2. Impact of various population effects on the M/L-age 
relation 

The effects considered in Sect. [4] and [5] which influence the 
integrated magnitudes and colours must also affect the mass- 
luminosity relation. As we discussed earlier the dynamical mass 
loss directly diminishes cluster mass. The scarcity of the popula- 
tion of stars in the upper Hertzsprung-Russel Diagram produces 
luminosity flashes of clusters. When one varies the lower mass 
limit of stars in a cluster with a given mass, the massive part of 
the IMF is changed. In this section we consider the manifestation 
of all these effects alone and in various combinations. 

In the left panels of Fig. [TT] we show models of different 
masses (M c = 10 3 , 5 • 10 3 , 10 4 m Q ) constructed with a con- 
tinuous IMF. Since the mass-to-luminosity ratio in the contin- 
uous case does not depend on cluster mass the different curves 
would overlap in the plot. To avoid crowding we shift (as indi- 
cated by thick arrows) the least massive model down, and the 
most massive model up from their original position. As one can 
see, the mass loss affects the later stages of the cluster evolution. 
At the same time clusters of different mass appear to have a dif- 
ferent state of dynamical evolution. For example, at t = 1 Gyr 
the model of M c = 10 3 itiq is the most dynamically evolved, 
and has considerably changed in the M / Ly compared to the 
model without dynamical mass evaporation. At the same time 
more massive models are practically insensitive to the mass loss. 
Note the good agree ment of our data with recent calculations of 
lAnders et al.l (|2009) (see their Fig. 5). In spite of a considerable 
impact on models of smaller cluster masses, the mass loss ef- 
fect is not able to change the evolutionary path of the cluster 
sufficiently: independent of the initial cluster mass, the M/L- 
age relation rises almost monotonically during the whole life of 
the cluster. Qualitatively the NIR-models demonstrate a similar 
behaviour as the optical ones. 

In the middle panels of Fig. fTTJwe show models of M c — 
10 3 m Q with various lower limits of stellar masses mi built with 
a continuous IMF. They also include the effect of mass loss 
due to dynamical evolution. Unlike the left and right panels, 
the curves are at their true position without any shift. The value 
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Figure 11. The evolution of the mass-to-luminosity ratio in the optical (upper row) and NIR (lower row) and various population 
effects. The curves shown in the left and middle panels present SSP models with continuous IMF, in the right panels we compare 
continuous-IMF and discrete-IMF cases. The left panels show the influence of cluster mass loss due to dynamical evolution for 
clusters of different initial masses with stellar low-mass limit mi = 0.1 m Q . The curves are artificially shifted vertically in order 
to avoid crowding as it is shown by cyan arrows. The upper (shifted up) curve corresponds to the initial mass of M c = 10 4 m Q , 
the middle (shown in the original position) curve is for M c — 5 • 10 3 to q , and the lower (shifted down) curve corresponds to 
M c = 10 3 m0. The dotted portions (as everywhere in this figure) mark models without mass loss. The middle panels show the 
effect of variations of the lower mass limit of stars for a model of M c — 10 3 ro . The upper, middle, and bottom curves correspond 
to mi = 0.01, 0.1, 1 rn©, respectively. In the right panels we compare models with continuous IMF (curves) and discrete IMF 
(connected symbols) for clusters of different masses, having mi = 0.1. Crosses indicate model with M c = lO 4 ™©, triangles 



M c = 5 • 10 3 m Q , and diamonds M c = 10" 
order to avoid crowding. 



mQ. As in the left panel, the curves for different M c are artificially shifted vertically in 



of mi varies from 0.01 m© (upper curve) through 0.1 m for 
the middle curve, and 1 rag for the lower curve. Again, dotted 
curves correspond to no-mass-loss models, and solid curves rep- 
resent evaporating models. The observed effect can be simply 
explained. In case of a model with mi = 0.01 m Q , the bulk of 
cluster mass comes from low-mass stars, making the population 
of luminous stars more rare implying a lower integrated lumi- 
nosity compared to the reference model of mi =0.1 m Q . This 
lead to an upward shift of the 0.01 m Q curve with respect to the 
reference model. In contrast, in a model with a lower limit of 
1 m no low-luminous stars will be formed, and the ratio will 
be shifted down with respect to the reference relation. In addi- 
tion, dynamical evolution will remove the bulk of cluster masses 
by evaporation of low mass stars. This process is more effective 
in clusters with low m;, and it will have almost no impact on the 
high-TO; counterpart. 

In the right panels of Fig. QT|we compare models of differ- 
ent cluster masses (M c = 10 3 , 5 • 10 3 , 10 4 m o ) constructed 
with a continuous and discrete IMF. In order to avoid crowding, 
we have again shifted artificially the curves of different masses. 
Here we observe a picture known from the previous sections: the 



discreteness of the IMF causes a spread around the continuous 
solutions. Whereas e.g. in B we practically do not observe any 
influence, a bias becomes noticeable in V, increases consider- 
ably in the NIR, and for K s the bias between the continuous and 
discrete models reaches a considerable value of about one order 
of magnitude for clusters with M c — 10 3 m . We also conclude 
that stochastic fluctuations in the M/L ratio increase when the 
cluster mass decreases. 



6.3. Comparison with observations 

In Fig. [12] we combine the plots shown in Fig [10] and Fig QT| 
We plot only the subset of models that gave the best fit to the 
data in the optical and NIR. The curves show the models for 
the cluster mass M c = 10 3 m Q with different low-mass lim- 
its (mi = 0.01, 0.1 and Itoq). The effect of mass loss due to 
the dynamical evolution is included in the models. By varying 
mi, the continuous models can partly cover the region occu- 
pied by observations in the optical. An additional scatter is in- 
troduced, especially for younger clusters, if one considers mod- 
els with a discrete IMF. However, the main reason for a large 
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Figure 12. Comparison of observed and theo- 
retical Al c /L, log i-relations for the cases of 
continuous-IMF (left panels) and discrete-IMF 
models (right panels), cluster mass of M c — 
10 3 mo, and variable lower limit of stellar 
mass mi. The upper row corresponds to the op- 
tical, the lower one to the NIR passbands. The 
designations for the empirical data are the same 
as in Fig. \10\ For the models dash-dotted blue 
lines (diamonds in the right panels) correspond 
to m; = 0.01 m.0, solid green lines (triangles) 
to mi — 0.1 mQ, and dashed red line (crosses) 
to mi = 1 rriQ. 



spread of the observed MjLy relations in Fig. [T2] is the un- 
certainty of individual data points (see Sect. 16.11 ). The typical 
evolution of the M/L ratio in the optical can, therefore, be bet- 
ter understood if we consider the evolution of the median values 
in Fig. [T2] For the youngest clusters, the best fit is achieved by 
evolutionary tracks with extremely low to; -limit. The clusters 
with 7.25 < logi < 8.25 are well described by the model with 
mi = 0.1 to , whereas the data for older clusters agree with the 
track with to; between 0.1 to q and 1 to q . We interpret this as a 
consequence of the mass loss in open clusters: the young clus- 
ters are formed with stars filling the full range of stellar masses, 
then, on a time scale comparable to the relaxation time, they start 
to evaporate (losing preferentially low-mass stars), and at older 
ages their effective low-mass limit sufficiently increases. This 
process appears to be even more dramatic than that predicted by 
the used N-body simulations we used to describe the dynami- 
cal evolution of open clusters. Probably, events like encounters 
with molecular clouds impact the stellar content of older open 
clusters in a very effective way. 

It is evident from Fig. [T2l (upper-left panel) that the contin- 
uous models are not able to fit the observations in the NIR: the 
models predict ratios too low compared to the observations, es- 
pecially for younger clusters. The situation is improved if one 
considers models with a discrete IMF in the upper-right panel of 
Fig-El In this case, one can recognise a tendency similar to that 
we observe for the optical passband. Young clusters tend to be 
fitted better by models with low m; -limit, whereas evolved clus- 
ters are described more adequately by the tracks with high m;- 
values. In contrast to the optical band, this picture is disturbed 
by strong fluctuations of the Mj L ratio towards low values due 
to the discreteness of the IMF. Their contribution to the observed 
spread of data points is comparable with the scatter introduced 
by randomly distributed errors of observation. 

We conclude that the evolution of the observed mass-to- 
luminosity ratio indicates a mass loss in open clusters with in- 
creasing age. In the optical, the large dispersion of M/L ratio 



is caused mainly by uncertainties of the data, though the IMF- 
discreteness favours the observed scatter (the effect increases 
with decreasing cluster mass). In the NIR, the IMF-discreteness 
is an important contributor to the Mj L dispersion. The higher 
sensitivity of the NIR fluxes to the number of red giants in clus- 
ters explains why the NIR-diagrams show a decreasing trend 
with cluster age. 



7. The flashing population of local clusters 

As we show throughout this paper, the discreteness of the stel- 
lar IMF plays an important role in the description of the inte- 
grated photometric properties of open star clusters and defines 
the flashing activity of the cluster population. In this section we 
try to formalise the "discreteness term" by deriving a condition 
for an emerging red flash in a star cluster. On the basis of this 
condition we also illustrate how the flashing activity of the clus- 
ter population depends on cluster parameters such as M c and 

TO/. 

From Fig. [9] (left panel) we conclude that "RG events" are 
prominent in clusters younger than logi w 8.6 where the MS- 
branch of the colour-age relation can be clearly recognised. 
Beyond that age, the MS-branch changes its slope and trans- 
forms to the RG-branch i.e. contains mainly clusters with red 
giants among the brightest stars. Therefore, in the following dis- 
cussion we consider only relatively massive stars, which have 
a chance to become a red (super)giant during the first 400 Myr 
(logi < 8.6), i.e. have masses above ~ 3to . 

Below we construct the condition for the appearance of a 
red flash event: when a cluster becomes considerably brighter 
and redder than in the Main Sequence stage, and after a certain 
rather short red phase returns back to the previous stage. This 
event is related to the transformation of the most massive star of 
a cluster to a red giant, and can be described in terms of stellar 
evolution. 



A.E. Piskunov et al.: The evolution of luminosity, colour and the mass-to-luminosity ratio of Galactic open clusters: 



13 



For the red giant flash to be noticeable, the total (MS + 
post-MS) lifetime of the most massive star in a cluster must be 
shorter than the Main Sequence lifetime of the next massive star. 
Otherwise, the second star simply joins the first one in the Red 
Giant domain, and the cluster remains red. This condition allows 
a cluster after a post-MS flash to become a MS-cluster again, as 
long as the second star still stays on the MS. In other words, if 
the most massive star has a mass to + dm, and the second mas- 
sive star has a mass m this condition can be written as 

TMs(m) > T ( m + dm) , 

where r and tms refers to the total lifetime and the main se- 
quence lifetime of a star of given mass, respectively, and dm <C 
m. The total and Main Sequence lifetimes differ by the life- 
time in the Red Giant stage r = r^/g + trg ( we assume that 
the RG-stage accumulates all the later phases of stellar life), or 
t ms = (1 — A) r, wh ere A = trq/t. As evolutionary cal- 
culations show (see e.g. [Bertelli et al. 1994: iGirardi et al] l2000: 
iMarigo & Girardil 120071: iMarigo et al.ll2008l) . A is a weak func- 
tion of stellar mass with an approximate value of A « 0.3 ... 0.1 
in the mass range log m/m Q = ... 1.9. Expanding the right- 
hand side of the inequality above into a Taylor series, one gets 

(1 — A) t(to) > r(m) + r (to.) dm + . . . 

Following stellar evolutionary calculations (e.g. Padova mod- 
els), the lifetime can be expressed as 

t(to) = /3 771~ 7 , 

where 7 > varies in the mass range log(m/TO Q ) = . . . 1.9 
from 3.8 to 0.3. Since dm is relatively small, we can assume, for 
simplicity, the parameters /3, 7 and A to be constant. Then the 
derivative can be computed as 



t (to) = — t(to) 



7 



and considering the linear terms only, we express the inequality 

as 

dm A 

— > - , (3) 



or in other terms 



7 



d log m > — log e . 

7 



This supports our previous conclusions in Sect.5 that the more 
massive the cluster is, i.e. the denser the cluster isochrone is pop- 
ulated with stars (d log m being small) the harder it is for the 
cluster to produce the flashing. This means, also, that the ability 
of a cluster to demonstrate red flashing is controlled by stellar 
evolution parameters: the duration of the red-giant phase A, and 
the decrement of the lifetime 7. The more pairs of stars having 
masses log to/to© > 0.5 will fulfil the condition (O at birth, 
the more RG-flash events a cluster will undergo in its life. We 
use this property to build a simple model of cluster population 
below. 

We stress that the condition $3^ does predict only the case 
of the colour bias with respect to the continuous SSP and does 
not deal with colour fluctuations around the SSP that can oc- 
cur even if (O does not hold. Thus, we consider this condition 
as a colour-bias constraint and apply it below to reproduce the 
flashing statistics of the local cluster population. We apply the 
colour-bias condition to a mass spectrum of a newly formed 
cluster constructed by a numerical simulation as described in 
Sect. 13. 31 The basic parameters of the simulation are cluster mass 



A 




- 










log M c 

Figure 13. The dependence of frequency n(RG)) of red flashes 
on cluster parameters as derived from Monte-Carlo experiments. 
The left panel shows a dependence on M c with 77i/=0.01, 0.03, 
0.06, 0.1, 0.25, 0.5,1 to q (from top to bottom along the dotted 
line). The right panel shows a dependence of n(RG) on to/ with 
M c =5 ■ 10 2 , 10 3 , 2.5 • 10 3 ,5 • 10 3 , 10 4 , 2 • 10 4 , 5 • 10 4 to q (from 
top to bottom). The bars indicate statistical uncertainties derived 
in the experiments. The horizontal solid line marks the observed 
frequency of local RG-clusters within the completeness area of 
the cluster sample of 850 pc. The hatched area corresponds to the 
Poisson error of this estimate. The vertical dashed line in the left 
panel shows th e estimated value of t he average mass of formed 
clusters from Pis kunov et alj d2008al) . 



M c , the lower limit of stellar masses mi, and the slope a of the 
IMF (assumed to be 2.35 across our study). These model pa- 
rameters define the population of the cluster and in particular, at 
masses greater than a limit of log too / to© =0.5 adopted above. 
For each pair of stars in this mass range, we compute the differ- 
ence in log to and apply the condition (O to them. Those stars 
which satisfy this condition are marked as flash-producing. In 
order to get statistically reasonable results, we repeat this proce- 
dure several times (usually a few tens of simulations are suffi- 
cient to get statistically meaningful results). As we noted above, 
we do not need to follow the whole evolution of a cluster to get 
the flashing-activity history. It is sufficient to analyse the initial 
mass spectrum only. 

The results of the simulations are presented in Fig.[T3lwhich 
shows the RG-flash frequency (the ratio of the number of stars 
satisfying condition (0 to the total number of stars formed with 
to > too) for different M c and to/. At M c > 500 to q the curves 
show a regular pattern: the RG-flash frequency decreases with 
increasing M c . At M c = 5 • 10 4 m Q , the cluster isochrone is 
so densely populated with massive stars that the RG-flashes are 
practically no longer possible at either value of to/. In contrast, 
at cluster masses smaller than M c = 500 to©, the flashing activ- 
ity of clusters is constantly high (n(RG) ~ 50 . . . 60%), but the 
isochrone is so poorly populated that the behaviour of the model 
curves is rather volatile. In the regular part of the diagram, one 
sees a diminution of the frequency n(RG) at given M c with in- 
creasing to/. Again, this occurs due to a rise of the stellar density 
of stars with log too / toq > 0.5 in the corresponding range of 
the isochrone. 

We also show in Fig. [T3] the empirical estimate of the flash- 
ing frequency in the Solar Neighbourhood which is found to be 
25 ± 4%. The estimate is derived from our ASCC-2.5 sample 
of open clusters which are younger than log/; = 8.6. To avoid 
a selection effect, we involve o nly clusters res i ding w ithin 850 
pc from the Sun. According to Piskunov et al. (2006), this dis- 
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tance corresponds to the completeness distance of our sample, 
and guarantees that at smaller distances we are free from the in- 
completeness bias. The estimate was computed as the ratio of 
the number of clusters having red giants (45) to the total number 
of selected clusters (183). The empirical value of n(RG) can 
be somewhat overestimated, since we are not sure that all ob- 
served RG-clusters undergo currently a flash (i.e. they may come 
back to the MS state after the observed RGs will die). However, 
as a careful inspection shows, most of the nearby clusters are 
"flashing-biased", i.e. are poor enough and contain only a few 
red giants. In the right panel of Fig.[13]we show the n(RG)—mi- 
relation with M c taken as a parameter. One can see that the ob- 
servations are compatible with initial cluster masses in the range 
M c = 10 3 ...5-10 3 to q . 

In principle, Fig. [T3l shows how statistics on the red giants 
(the brightest stars of a cluster) can be converted into a knowl- 
edge on the lower mass limit of forming stars provided that the 
initial cluster mass is known. The typical mass of young clus- 
ters of 4.5 • 10 3 TTtfD we found from t he analysis of the cluster 
mass function ( Pis kunov et al . 2008a) indicates that to/ should 
be small. Taking into account the uncertainty both in the fre- 
quency, and in the initial mass values, we expect a possible to/ 
range of 0.01 . . . 0.1 m . 

8. Summary and Conclusions 

We have proved the necessity to modify the traditional SSP- 
model widely applied in studies of the extragalactic clusters be- 
cause it fails to explain the luminosity and colo ur distribution 
and evolution of Galactic open clusters (Paper II). The observed 
colour offset between theory and observation led us to the idea to 
revise the basic assumption used in present-day SSP - the conti- 
nuity of the IMF of clus ter stars . The successful solution of the 
observed discrepancy in lPaper 11 provided an impetus for the ex- 
tension of the investigation of the properties of the discrete SSP 
model. This paper includes the analysis of the evolution of lumi- 
nosity, colour, and the mass-to-luminosity ratio with cluster age 
of both continuous and discrete SSP-models. The study includes 
several effects which either occur only in a discrete approxima- 
tion (like the mi -effect in cluster colours), or were rarely studied 
even in the continuous approximation (mass loss from clusters). 

The evolution of the integrated magnitudes shows that mass 
loss is the weakest effect among all others taken into account. 
It caused only a small change of cluster magnitude of older 
clusters. This can be expected since in the considered regime, 
mass loss occurs mainly due to evaporation of the low-mass stars 
that only weakly contribute to the cluster luminosity. However, 
the mass loss effect taken into account modifies the general be- 
haviour of the evolutionary track, increases the steepness of the 
magnitude-age relation, and thus provides a better agreement 
with the observations. The effect can possibly become more 
prominent when the other causes of mass loss (which are not 
considered here) are taken into account, e.g. encounters of clus- 
ters with molecular clouds. 

The effect of the lower limit of stellar masses to/ is stronger. 
Being applied to a cluster of a particular mass it alone is able 
to account for the total width of the observed magnitude-age re- 
lation. Unlike the mass loss it does not affect the slope of the 
relation. 

Last but not least, the IMF-discreteness effect has large con- 
sequences for the cluster luminosity. It is responsible for cluster 
flashing and track fluctuations due to the appearance of red gi- 
ants. Moreover, it is the reason for a systematic bias between lu- 
minosities of MS-clusters and the predictions by standard SSP- 



models. The effect is especially strong at long wavelengths, 
where the flash amplitude can reach 4 mag, and the branch of 
MS-clusters may be located 2-3 magnitudes below the continu- 
ous track. 

In summary, we conclude that the general evolution of 
Galactic open clusters in the magnitude-age diagram is ade- 
quately described by models of moderate mass (5 ■ 10 2 ... a few 
of 10 3 m Q ). We explain the difference between th is value and the 
averag e mass of forming clusters of 4.5 • 10 3 m Q (Piskunov et al. 
2008a) by the episodes of cluster mass loss (We idner et al. 
120071) . which occurs in young clusters as a result of violent 
events (SN explosions, HII formation, and stellar winds), by 
mass loss due to encounters with molecular clouds, and due to 
regular evaporation. 

All this indicates that the approximation of the observed re- 
lations with a simple grid of continuous SSP-tracks is much too 
simplified. For example, the extragalactic cluster samples can be 
well contaminated at magnitudes My = —6 • • • — 10 mag with 
flashing clusters, which in reality have lower MS-luminosity and 
cluster mass. This in turn can make the observed sample artifi- 
cially overpopulated by luminous and massive clusters and will 
make the apparent luminosity and mass functions flatter. 

Therefore, the direct implementation of standard SSP mod- 
els, as it is usually applied to extragalactic clusters to estimate 
their masses, must be taken with caution. The only possibility 
to avoid large errors in mass estimates is using observations at 
short wavelengths (e.g., U, B). 

Cluster colours also stress (sometimes even more than the 
magnitudes do) the trends mentioned above. For example, the 
discrete models demonstrate a colour-dependence on the low- 
mass limit, or on M c . Sometimes we see that both effects mimic 
each other in the discrete case. This is impossible in the case of 
models with a continuous IMF where the cluster colours are, by 
definition, independent of M c and to/ . 

Again the basic difference between the standard and discrete 
models is the colour bias produced by MS-clusters, and the flash- 
ing activity related to the formation of red (super)giants. As for 
integrated magnitudes, the prominence of both events depends 
on the population of the massive end of the stellar mass func- 
ti on in a cl uster and is directly related to M c and to/. As shown 
in iPaper ij, with an average value in A(B — V) 0.3 mag 
the colour bias reaches its maximum of the order of 0.5 mag at 
log i«7 and can be still detected at log t > 8. Again, the colour 
bias increases at longer wavelengths and reaches a value of the 
order of 1 mag in (J — K s ). The presence of the bias in the stan- 
dard SSP-models should be taken into account by determining 
cluster reddening and age from the integrated colours of clusters 
with masses less than 10 6 m^, and warns against a straightfor- 
ward use of the SSP-grids (IPaper ij) . The discrete model with a 
cluster mass of 10 3 ... is successful in reproducing the colours 
of young clusters (\ogt < 8.5) of our cluster sample. At older 
ages the flashing population is represented mainly by clusters 
of lower masses with different evolutionary stage of red giants 
along the RG-sequence. 

It is evident that an independent approach for the determi- 
nation of reddening and mass is required before cluster colours 
can be used for the age estimate. But even in this case, the dis- 
creteness effect must be taken into account, and the age will be 
determined within some confidence interval by interpolating a 
grid of a properly selected SSP-model. In order to increase the 
reliability of the results, the age determinations should be based 
on observations in short wavelengths. 

We have constructed the mass-to-luminosity ratio based on 
luminosity and mass data determined independent of each other 
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for Galactic open clusters. This makes the derived relation M /L- 
age bias free. The average relation increases with age in the op- 
tical, and decreases in the infrared. The observed scatter turns 
out to be slightly larger than that expected from the accuracy 
of the input values in the optical, but considerably larger in the 
infrared. 

If we introduce the effects described above in our models, 
the following picture emerges. Mass-loss in clusters makes the 
M/L tracks to appear flatter at later stages (log/; > 8.5). The 
less massive the cluster is, the more pronounced is the effect. The 
combination of mass loss and mi makes the curves almost flat in 
the NIR after log t m 8. The variations of mi within the limits 
adopted in this study are able to produce a considerable vertical 
scattering in the relation, and the IMF-discreteness adds extra 
fluctuations. In the NIR at smaller cluster masses one observes 
a strong bias in the relation, implying that the continuous option 
predicts too low M/L. 

Compared to the observations, the tracks show better agree- 
ment with a discrete IMF, covering practically the complete ob- 
served M/L spread. One also notes an interesting detail, the 
younger (log/ < 7.5) clusters are in better agreement with 
tracks with lower values of mi, and older clusters (log t > 8.5) 
agree with higher-m; tracks. We interpret this finding as an indi- 
rect evidence of the dynamical evolution of open clusters, born 
with a full mass spectrum of stars, which is narrowed from the 
massive end by stellar evolution, and from the low-mass end by 
dynamical evolution. Probably, the mass loss due to irregular 
events (e.g., encounter with molecular clouds) may play an im- 
portant role on the decreasing of the M/ L-ratio in older clusters. 

Taking into account the time-scales of stellar evolution, we 
derived a condition for the occurrence of colour and magni- 
tude flashes observed in the photometric diagrams for open 
clusters. This condition reflects well the changing properties 
of stars caused by t heir e volution. Unlike the approach by 
ICervino & Luridiana (2004) to separate physically feasible and 
un-realistic SSP-models, our condition reveals the origin of a 
bias in realistic models. Our model appears to be very convenient 
in predicting the flashing activity of clusters. We use this "flash- 
ing condition" for the construction of a simple population model 
of star clusters, which relates the frequency of RG-flashes, clus- 
ter mass, and low-mass limit of cluster stars. The application of 
this model to open clusters in the Solar Neighbourhood has the 
following consequences: with the average observed frequency of 
flashing clusters of 25 ±4% and with an average mass at birth 
of a few 10 3 m0, we infer that in newly formed clusters the 
lower limit of the stellar mass spectrum should fall into the range 

0.01... 0.1 7710. 
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